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Abstract. We discuss accretion flows on massive black holes in which different elements of the flow (clumps) have velocities 
that may differ substantially. We estimate the consequence of collisions between these clumps as they come close to the central 
object and calculate the resulting radiation. We show that this radiation is similar to that observed in the optical to X-ray spectral 
domain in Seyfert galaxies and quasars. We also show that the large scale accretion is likely to be clumpy when arriving in the 
active region and that the clumps keep their identity between collisions. 
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1. Introduction 

Although the basic elements of the physics of Active Galactic 
Nuclei (AGNs) were established soon after the discovery of 
quasars in the early 1960's and are now beyond dispute, 
many aspects of their physics still have to be discovered. 
We know that AGNs are powered by the gravitational en- 
ergy of matter that falls into the deep gravitational well cre- 
ated by massive black holes, and can estimate the mass of 
the black holes and the accretion rates (or at least lower 
limits to these quantities) from the bolometric luminosity of 
the objects. The masses of the black holes may also be de- 
duced from kinematic measurements (e.g. [Kaspi et al. 2000| 
IPeterson et al. 2004llPaltani &Turler 2005> . We have also iden- 
tified many of the emission processes that create the continuum 
spectral energy distributions and the emission lines. We can 
deduce the state of some of the emission gas from line ratios 
and the shape of the continuum emission. These facts are now 
well documented in a number of books and reviews (see e.g. 
ISetti & Swings 2001> . 

While our understanding of the radiative processes that cool 
the electrons is very advanced, our understanding of the pro- 
cesses that transfer energy from the accreting mass to the ra- 
diating electrons is still quite unsatisfactory. In a stationary 
situation, the electron cooling rate must be equal to the elec- 
tron heating rate. There are several ways to transfer energy 
to electrons. Electric fields do so. It is difficult, however, to 
see how organized electric fields could arise across 100 Rs 
(Rg = M being the mass of the black hole) from the 
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energy of the accreting ions. Coulomb collisions are also a pro- 
cess in which energy from ions is transferred to electrons. We 
will see in Sect.lsjone case where this process plays an impor- 
tant role. 

Optically thick and geometrically thin accretion disks 
( |Shakura & Sunyaev 1973|/ have been used widely either to 
understand accretion and the corresponding radiation of the 
gravitational energy in the blue bump of AGNs or as sup- 
port of more complex structures including Comptonising coro- 
nae (Haardt et al 1994i. This continues despite the difficul- 
ties met by accretion disk models in accounting for the vari- 
ability in the blue bump of AGN (see Sect. |2ji and the lack 
of accepted sources of viscosity in the disk material large 
enough to remove angular momentum at the needed rate 
JKoratkar & Blaes 1999> . 

In the following we propose a framework in which the 
accreting flow is much more structured and clumpier than a 
standard accretion disk, and we then make explicit use of the 
dumpiness of the flow to explain the observed properties of 
AGN. We consider a flow in which shocks are the primary 
means through which the accreting mass transfers its kinetic 
energy to the radiating electrons. We find that cascades of 
shocks are expected, first in an optically thick environment and 
then in optically thin conditions. We estimate the properties of 
the radiation emitted and show that they match those observed 
in the optical to X-ray emission of AGNs. We can also un- 
derstand the origin of relativistic electrons in this framework. 
However, the presence of very organized relativistic outflows 
(jets) does not follow from the shocks we discuss, but requires 
an additional phenomenology, probably based on the rotation 
of the central black hole. 
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We consider accretion flows characterised by very inhomo- 
geneous conditions and use the fact that shocks in such a flow 
are an efficient way to transform bulk kinetic energy in thermal 
and non-thermal relativistic particles. We then calculate the re- 
sulting radiation field. 

In Sect. |2] we summarize the most prominent observational 
evidence on AGN that are relevant for our discussion. We then 
present the main features of the model and the nature of the 
accretion flow in Sect.|3] In Sect. |4]and Sect. |5lwe describe 
quantitatively the shocks in the accretion flow that energetize 
the electrons and discuss the resulting electromagnetic radia- 
tion. 



2. Some observational results 

We summarize here some observational results, mainly based 
on 3C 273, emphasizing those aspects which are difficult to un- 
derstand in the frame of accretion disk models with or without 
a corona. 

2. 1. Blue bump 

Seyfert galaxies and quasars have a broad optical-UV emis- 
sion that is markedly in excess of the extrapolation of 
the infrared emission to high frequencies. This compo- 
nent, the blue bump, has been associated with the accret- 
ing matter since Shields (1978) This emission is broader 
than a single black body, and has been modeled as an 
optically thick and geometrically thin accretion disk by 
jMalkan & Sargen t (1982)| and many authors since then. This 
component has been associated with shocks in an outflowing 
wind by' Camenzind & Courvoisier (1983) ] Despite many years 
of efforts, the blue bump has still not been described well by a 
standard optically thick and geometrically thin accretion disk, 
even when several effects that may distort the original disk 
emission are taken into account (Koratkar & Blaes 1999 1. 

The shape of the blue bump, and in particular the energy 
at which it is observed, is very similar in objects of very dif- 
ferent luminosities ( Walter et al. 1994 ). Accretion disk models 
predict, however, that the temperature of the disc depends on 
the mass of the object, hence on its luminosity. The lack of cor- 
relation between the shape of the blue bump emission and the 
object luminosity can thus be compatible only with standard 
accretion models if the mass of the black hole and the rate of 
accretion are in a finely tuned relation ( .Walter & Fink 1993 ) . 

It has been known for a long time 
( ICourvoisier & Clavel 19911 and references therein) that 
the blue bump varies so that there is at most a very short lag 
between the UV and the optical emission. This suggests that 
the information about the flux variations travels at speeds close 
to the speed of light between the hotter and the cooler parts of 
the blue bump emission region rather than at sound or viscous 
velocities. 

The optical and ultraviolet variability of 3C 273 has been 
studied by iPaltani et al. (1998)| It was found that the lag be- 



tween the 1250 A and long wavelength lightcurves. At (A), de- 
pends on the wavelength in the following way; 

At(A) ~ 4days • log2(A/1275A). (1) 

A study of N GC 7469 by [Wanders etal.(1997)| and 
Colli er et al. (19 98) also showed a delay increasing with wave- 
length and suggested that the delay in this object is about 
0.6 days rather than 4 days as in 3C 273. 

It has been shown by Paltani et al. (1998) that the UV Ught 
curve of 3C 273 can be described by the superposition of in- 
dependent finite duration events. The available data constrain 
the shape of the individual events only crudely, which can be 
characterized by a linear increase up to a maximum luminosity 
followed by a linear decrease. The increase and decrease du- 
rations may be assumed to be the same. This parametrization 
is also used in the following. It will be shown that the details 
of the event shape do not influence the conclusions that can be 
reached. The rate of events was found to be about 3 per year, 
their duration 2/i (fi being the time needed to reach peak inten- 
sity) to be about 0.4 year, and the energy of each event to be on 
the order of 10^^ erg. 

2.2. X-ray - blue bump correlations 

The X-ray slope in the 2-10 keV domain of 3C 273 is about 
0.5 and shows some variations. INTEGRAL observations in 
January 2003 showed that at this epoch the hard X-ray emis- 
sion (2-200 keV) was at one of the weakest historical levels and 
that the slope was considerably steeper than at earlier epochs 
(ICourvoisier et al. 2003> . 

The short lags between the light curves of the blue bump 
emission at different wavelengths has most often been ex- 
plained by assuming that the gravitational energy released by 
the accreting matter is not dissipated in the optically thick ac- 
cretion disk, but instead in a corona outside the disk that illumi- 
nates it and heats it. In this way, variations in the heating (hard 
X-ray) flux indeed causes quasi-simultaneous variations in the 
different parts of the blue bump emission material. This model 
predicts that the heating flux and the reprocessed blue bump 
emission are tightly correlated at zero lag. Measurements of 
the heating flux are difficult, as most of the heating energy is 
radiated at energies higher than 10 keV, whereas most of the 
X-ray photons are emitted at energies less than lOkeV. Most 
of the observations were therefore performed at energies less 
than lOkeV. Observations of the quasar 3C 273 in the hard X- 
rays were performed with the BATSE instrument on CGRO and 
with INTEGRAL, they do not show the expected correlation 
with the blue bump emission(Fig. 1). 

Studying the cross correlation between the X-ray and the 
UV Ughtcurves of 3C 273, |Paltani et al. (1998b), observed that 
the cross correlation has two relative maxima, one at zero lag 
and the other suggesting that the UV emission leads the X-ray 
emission by approximately 1.8 years. This result is based on 
observation of one major event that was well covered by both 
light curves, it will therefore need confirmation when many 
such events can be studied simultaneously. 
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Fig. 1. Comparison of a) the optical 
B band lightcurve and b) the 60 keV 
X-ray lightcurve of 3C 273. Optical 
data are from the Geneva photomet- 
ric system (open points) and scaled 
observations from Kaspi et al. (2000) 
(filled points). X-ray observations are 
from CGRO/BATSE (averaged into 30- 
days bins) (open points) and from 
INTEGRAL/ISGRI (filled points). 



2.3. Line emission 

Quasars and other AGNs (with the notable exception of BL Lac 
objects) show prominent line emission. The line emitting mate- 
rial has a density of about 10^" cm~^, and is photo-ionized (see 
e.g. JVetzer 1990) . This line emission material fills only a small 
fraction of the volume, it has been discussed mainly in terms 
of isolated clouds. The clouds have not yet been observed, nor 
has the discrete nature of the line emission been established. 
IDietrich et al. (1999)| attempted to measure the signature that 
the presence of discrete clouds should leave on the line pro- 
file in 3C 273. Even at very high resolution and high signal- 
to-noise, the line profiles are smooth. Only lower limits to the 
number of clouds contributing to the lines could be found. The 
number of individual clouds in the broad line region of 3C 273 
is thus larger than ~ 10*. This very high number may imply 
that a description in terms of individual clouds may not be ap- 
propriate, but that a better model may be found in a highly dis- 
organized matter flow. 



3. Cascades of shocks 

In order to explain the complex properties of the emission of 
AGN outlined above, we discuss the emission expected from 
an accretion flow that is characterized by individual clumps 
rather than a smooth flow. The flow generates first optically 
thick shocks, then optically thin ones. These shocks are are the 
origin of the radiation. We consider a clumpy flow in which 



the clumps have velocities given by the gravitational field of 
the central black hole and a wide distribution of velocity di- 
rections. Clumps following orbits in the gravitational field of 
the central black hole will lose energy either through tidal pro- 
cesses in the vicinity of the black hole or when they interact 
with one another. Only the latter process can dissipate the bulk 
kinetic energy of the clumps and thus be at the origin of the 
large luminosities of AGNs. We therefore consider in the fol- 
lowing the shocks that form at the boundary of two colliding 
clumps. We then follow the shocked material as it interacts with 
surrounding material having also been shocked in similar con- 
ditions. We show that the radiation emitted as a consequence of 
these processes closely matches the observed optical to X-ray 
properties of AGNs. 

The phenomenology we describe here relies on the exis- 
tence of a very clumpy accretion flow. This is expected on at 
least two grounds (i) the interstellar matter from which accre- 
tion in AGNs is nurtured is very structured and certainly not 
smooth; (ii) self gravity in the accretion flow is likely to play 
an important role at distances larger than 1 pc from the black 
hole (see ' ColUn & Zahn 199^ and references therein). These 
authors discuss how stars of about 10 solar masses form in the 
1-lOOpc region of an accretion disk and how they migrate to 
the central regions. The mass range of stars that are being pro- 
duced does correspond to the clumps at the origin of the phe- 
nomenology described here. 

In past years [C ourvoisi er et al. (1996)| and 
|Torricelli-Campioni et al. (2000), have proposed and dis- 
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cussed a model in which the accreted gas clumps were 
assumed to be stars. Although it may well be that a fraction of 
the accreted matter is indeed in the form of stars, the approach 
taken here does not rely on this assumption. It is sufficient to 
request that the primary shocks be optically thick. The nature 
of the clumps can thus be very diverse, from stars (be they 
from the host galaxy as in our previous work or generated in 
the accretion flow), to dense elements of an accretion disk. 
The clumps may also form in the accreted interstellar medium 
as discussed in Sect. or even in outflowing winds. Some 
properties of the clouds will be discussed in Sect.0 For shocks 
to form, it is mandatory that clumps have velocities that differ 
either in direction or magnitude so that the relative velocities 
are the same order as the clump velocities. This suggests that 
clumps have a broad distribution of angular momentum, as 
expected if the accreted matter comes from regions that have 
been disturbed by either non-spherical potentials or from 
merging and/or interacting galaxies. The colliding clumps 
may then be thought of as belonging to two (or possibly 
a distribution of) flows, each characterized by its specific 
angular momentum and origin. We consider here shocks 
between "generic clumps". The results obtained will then be 
used to develop more realistic approaches to AGN including a 
discussion of distributions of clump properties like mass, size 
and angular momentum, and a discussion of the properties of 
the region surrounding the black hole like density and velocity 
distributions. 

4. Optically thick shocks 

The energy content in the outbursts of the blue bump 
of the bright quasar 3C 273 has been calculated by 
jPaltani et al. (1998)| to be ~ 10''^ ergs . This is the kinetic en- 
ergy of some 10^'^ g at 0.1c, the Kepler velocity at 100 Rs- 
Since not all the kinetic energy is radiated, we propose that 
typical shocks involve clumps on the order of several lO'^'^g. 
Such masses agree with the discussion in Sect.|3] 

We therefore consider shocks between clumps of mass 
M — M^^IQ'^^ g and size moving at WKcpicr in the gravita- 
tional field of the central black hole. A/33 will typically range 
from a few to a few 10s. Around 100 i?5, fKcpicr = l^o.i 0.1c 
with (3o.i of order 1. We assume the clumps to be sufficiently 
dense for the photon diffusion time t^is within the matter to be 



very long compared to the crossing time tc 



"Koplo 



We 



We assume in the following spherical expansion for the sake of 
simplicity. Detailed numerical simulations of stellar collisions 
(Freitag 2000 1, using smooth particle hydrodynamics methods 
and including the internal structure of the stars, confirm that 
high velocity collisions between stars lead to the disruption of 
the stars and to a high velocity expanding gas cloud. 

The expanding gas will radiate a fraction 77 of its energy 
E at its photosphere. The time it takes for this energy to be 
radiated can be estimated by comparing the photon diffusion 
time tdiff with the expansion time. The photon diffusion time is 



/ ,R^2 

tdiS — - (t) = CTt JT-o, 

c I c 



(2) 



where I ~ (ctt • n)^^ is the photon mean free path, ric the 
electron density, R the size of the expanding cloud, and ctt 
the Thomson cross section. Use of the Thomson cross section 
is sufficient as it is expected that the shocks will completely 
ionize the material. The expansion time is tcxp = jr~- 

The expansion time and diffusion times are equal; i.e. 
^cxp ~ tdis, at i? = i?max, which may be expressed as 

i?max = — - CTT —) ' = S.lA/g;^ 10 ''cm, (3) 

c mp 47r 

where now /?o.i — Wcxp/ (0.1 c). i?„iax depends on the expan- 
sion velocity and the mass of the clumps. This size is reached 
after a time tcxp(^max) = , which for two mass clumps 
at 100 Rs, i.e. with WKcpicr — 0.1 c is on the order of 



<cxp = 1.0A4/'/3o:1^' lO^s 



(4) 



will consider optically thin shocks in sect.|5l 

When two such clumps collide, their kinetic energy in the 
center of mass frame is thermalized. The resulting configura- 
tion is a very hot dense gas cloud in adiabatic expansion. This T{t) = [ 
is similar to a supernova in many ways. The properties of the 
expansion are therefore expected to resemble those of expand- 
ing supernovae in their early stages. The total energy of the 
gas cloud in its rest frame, that of the center of mass of the 
two coUiding clumps, is E ^ M WK^pi^^j. f^goom, where ?7gcom 
is a factor taking the geometry of the collision into account and 
is assumed to be 1 in the following. Coulomb collisions at the 
origin of the thermalization being elastic, we expect that the ex- '^^^^ 
pansion velocity of the gas cloud is Wcxp — "J^Kopioi - The geom- „ , L,, 
etry of the expanding gas reflects the symmetry of the colhsion. 



This time is comparable to the typical rise time of an event 
/i ~ 0.2 years, as measured in the UV light curves of 3C 273 
(see Sect.B'Pa itani et al. (199S% for masses of several 10 e. 

The expanding matter is very optically thick. The luminos- 
ity of the expanding mass is therefore given by the size of the 
photosphere and its temperature, which depends on the rate of 
the heat transfer to the exterior. The calculation of this rate is 
beyond the scope of the present paper. We can, however, es- 
timate the temperature of the photosphere by assuming that it 
expands with Uoxp and that its luminosity is given by the sim- 
ple model of jPaltani et al. (1998)| in which the luminosity in- 
creases linearly until it reaches a maximum at a time t — n after 
the beginning of the outburst. We identify /i with tcxp(^^max)- 
With giving the luminosity at t = /i, one obtains for the 
time dependence of the temperature T{t), where t is the time 
elapsed since the collision between two clumps: 



11/4 



(5) 



where a is the Stefan-Boltzmann constant. 

The time dependence of the temperature (Eq. |5|l may be 
written as 



T = 



Y'^ = 1.1 X io4( 



J44 



cxp 



M33/30.: 



.)1/4k 



(6) 



(7) 
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is the temperature reached at the photosphere at i = /i and 
L4A = Yp^ergss^^. The event properties derived from the 
observations of 3C 273 by Palta ni et al. (199}^ suggest that the 
luminosity is a few 10"'^ erg s~^. 

This estimate is valid only provided that the temperature of 
the material in expansion is sufficient for the surface to emit at 
the blackbody limit. We can check that this is the case by es- 
timating the temperature at the begining of the expansion and 
the adiabatic rate of decrease in the temperature. The initial 
temperature is given by the shock velocity (~ 0.1c) and is 
of the order of 10^" K. The temperature decrease is given by 
the adiabatic expansion of a radiation pressure dominated gas 
T ~ Thus even assuming an initial clump or shock size 
as small as that of a star (10^^ cm), the gas expanding to i?niax 
does not adiabatically cool below the photosphere temperature 
estimated above. 

The model proposed here is able to account for the blue 
bump emission of AGNs very well, the typical temperatures 
of which are indeed on the order of 10"* K and dominates the 
optical UV emission. The model also predicts a broader com- 
ponent than a single black body, because several events will be 
co-existing in the nucleus at different stages of their evolution 
and therefore at different individual temperatures. This matches 
the observations. Different objects are expected to have blue 
bumps at roughly the same wavelengths, as this is given by 
collision parameters and not by the properties of an accretion 
disk that depends on the mass of the central object and the ac- 
cretion rate. Differences between objects will furthermore be 
small, because of the low power in which the event parameters 
enter Eq.|6l These parameters are such that the energy radiated 
during the event is on the order of 10^^ ergs. 



black holes. With an event rate of several per year, the gas en- 
velopes of the clouds therefore overlap in the central regions. A 
large fraction of the expanding gas of any cloud will be shocked 
again, this time, however, in an optically thin environment. The 
relative velocities can be expected to be somewhat larger than 
the sum of the Kepler velocities at the originating optically 
thick shocks due to the fact that gas moving inward and gravita- 
tionally accelerated has a higher probability of interacting with 
other expanding clouds than the gas moving outward. 

The temperatures reached at the optically thin shocks are 



kT 



16 



(9) 



Here, v^d is the relative speed of the expanding clouds where 
they shock and m is the average mass of a particle (=1 10^^"* g 
for completely ionized matter). With of about 0.1c, the 
temperatures obtained are on the order of 1 MeV. 

The temperature of the electrons is given by the equilibrium 
between their cooling and their heating rates in the shocked 
material. The electrons cool predominately through Compton 
cooling. They gain energy through Coulomb interactions with 
the protons and nuclei, which do not cool through Compton 
radiation (or at least do so much less than electrons). The 
Coulomb interactions of the electrons with the protons and nu- 
clei is a relatively slow process. One therefore expects that the 
electron temperature is less than that of the protons. 

The average electron energy can be estimated by consid- 
ering the equilibrium between energy gain through Coulomb 
collisions and Compton losses. Compton losses are given by 
(e.g. |Rybicki & Lightman 1979> : 



4. 1. Lags in the light curves 

Equation |6l can be inverted to estimate the time it takes for the 
temperature of an event to cool from Ti = T{ti) to T2 = 
T{t2). With AT = Ti - T2 one obtains the following for AT 
small compared to Ti 



J-i Ji 



(8) 



to the first order in AT/T. 

It therefore takes about 10*^ s to cool by 25% at T^. 
This time is very similar to the lag At (A) measured by 
IPaltani etal. (1998)lin 3C 273 between the visible and the UV 
light curves (see Eq.[Tll. Furthemore, the time delay increases as 
the temperature difference increases, corresponding to the ob- 
servation that the lag of lightcurves increases as the wavelength 
difference between the light curves increases (Eq.Q. 



5. Optically thin shocks 

The gas emerging from the expanding regions we have just de- 
scribed moves ballistically with bulk velocities Vcxp close to 
the Kepler velocity at the corresponding radius. Since the gas 
expands beyond i?max, it becomes optically thin at times larger 
than tcxp(^^max) and fills a volume of : 



that is com- 



parable to the region within lOORs, even for massive central 



iCompton = -^CTTCj'^ P'^ ' "ph, 



(10) 



Wph being the photon energy density and the other symbols hav- 
ing their usual meaning. The photon energy density is the one 
generated by the optically thick shocks of luminosity Lthick: 



ick 



AttRL 



29- 



L 



44 



M33 f3o. 



-erg cm 



(11) 



The electrons are heated through Coulomb collisions with 
the hot protons and ions at a rate Lcouiomb = Ep/tc where 



Ep is the proton kinetic energy (c 



1 MeV) and tc 



is given by the standard Coulomb energy transfer time (see e.g. 
Benz 2001 1: 



(12) 



4Ad • vl 



with 

STTneHnjA) ln(A) 

= 5 ==3.210 -ns-^— 

mj 20 

in c.g.s. units. The density may be estimated by 



-M™i^^5.3xl0^Af3-/^/34/Lcm-. 



47r/3i?: 



3 

max 



(13) 



(14) 
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Using the photon energy density from Eq. ([HJ and the 
electron density from Eq. il4i . the equihbrium icompton — 
Lcouiomb can be solved for the electron energy: 

~ Q.bMli' 13^\'^ L-J" Ep,M.Y. (15) 

where E'p^McV — i?p/lMeV. This estimate gives an average 
electron energy of some 250 ke V, which is less than the proton 
temperature, as expected. It depends with a low power on the 
parameters of the model. 

Shocks are well known for producing relativistic particles 
that contribute several percent of the energy balance (see e.g. 
IKirk 199 4* for a discussion of particle acceleration in shocks). 
This will also take place in the environment we have just de- 
scribed. In the presence of a magnetic field, possibly also ac- 
creted with the in-falling matter, these relativistic particles will 
produce synchrotron emission. 

5. 1. X-ray emission 

The Compton process leading to the cooling of the ~ 250 keV 
electrons will radiate X-rays. For a non relativistic electron 
population the spectrum is characterized by a power law of in- 
dex ( |Rybicki & Lightman 1979| l a = — , where r is the 
optical depth of the region. Using the density of Eq. (I14t and 
^max as the size, one can estimate that r ~ 3.6PqI. A is the 
relative energy gain per Compton diffusion, which is given by 
A = AkTc/rricC^ for non relativistic electrons. kTc is the av- 
erage energy of the electrons as estimated in Eq. ( I15> . For the 
parameters of our model, a ~ —0.3. This slope agrees with 
the observed hard X-ray slope of 3C 273 of about -0.5, con- 
sidering that the spherical geometry used to deduce the optical 
depth, and hence the spectral slope, is certainly a very rough 
estimate at this stage. 

It can be seen from Eq. il5\ that an increase in the blue 
bump photon energy density leads to a decrease in the electron 
temperature. This in term leads to a decrease in A and to a 
steepening of the X-ray emission. This behavior corresponds to 
the observations of 3C 273 made in 2003, when the blue bump 
was found to be at a very high level and the INTEGRAL slope 
was one of the steepest ever observed ( Courvoisier et al. 2003 1. 

We can estimate the X-ray luminosity of the optically thin 
shocks starting from the Compton luminosity icompton of a 
single non relativistic electron (Eq.fTOJ) with the photon energy 
density from Eq. and the electron velocity deduced from 
Eq. O: 

iCompton ^ 7.7 10-13m-6/7^-8/7^5A grgss-^. (16) 

For a shock in which there are /nip electrons the X-ray 
luminosity is 

Lx ^ 4.6 m^^Ml^'p-'/'Lli'cTgss-K (17) 

Time variations of the resulting X-ray emission are ex- 
pected on many time-scales that reflect the inhomogeneities of 
the colliding gas. These inhomogeneities result from those of 
the original gas as well as from those, on larger scales, due to 



the brightening and fading of the corresponding optically thick 
shocks. 

Correlations of the X-ray emission with blue bump light 
curves are expected at short lags, as the soft photons are be- 
ing Comptonised, and also at longer lags, because the opti- 
cally thin shocks at the origin of the hot electrons have evolved 
from the optically thick shocks that created UV emission. The 
lag in this second case is on the order of several months (> 
ioxp(^max)), the UV light curve leading the X-ray emission. 
This behavior corresponds well to the observations described 
in |Paltani et al. (1998b) where it is shown that the cross cor- 
relation of the X-ray and UV light curves of 3C 273 has two 
maxima, one indicating that the UV light curve leads the X- 
ray emission by ~ 1.75 year and the other indicating a good 
correlation with small lags. 

In the reference frame linked with the central black hole, 
the shocked material will have an angular momentum that is 
the vector sum of that of the material of the two original gas 
clumps. The shocks thus provide a means of reducing the angu- 
lar momentum per particle and thus will allow the shocked ma- 
terial to fall deeper in the potential well, where further shocks 
may be expected. These cascades serve as an effective viscosity 
in the accretion process, while they contribute to transforming 
the gravitational energy of the incoming nucleons to the radiat- 
ing electrons. 



6. Line emission 

The gas moving outward beyond Rmax, rather than inward fol- 
lowing the expansion of the optically thick shocked material 
described in Sect.|5j will be slowed down to velocities less than 
0.1 c by its ballistic motion. This material is expected to be very 
inhomogeneous and have a wide range of velocities. Its density 
is given by Eq. ( I14t . which corresponds to the densities ex- 
pected in the broad line emission region. We therefore suggest 
that the line emitting gas is the result of the shocked material 
once it has expanded in the outer regions of the nucleus. 



7. Ciouds 

The structure of the clouds and the rate at which they initially 
fall into the central region of the AGN are given by the condi- 
tions in the central regions of the host galaxy and the geometry 
of the gravitational field which drives the accretion from the 
large galactic scales to the immediate vicinity of the black hole 
where the high energy phenomenology is generated. 

In order to investigate how clumpy the accreted matter is 
likely to be, we consider the large scale accretion flow around 
an object accreting matter with a radial velocity Vacdf) = 
XV f f (r) at the Eddington rate, where Vff (r) is the free-fall ve- 
locity onto the black hole. Those clouds that form while falling 
into the central regions of the AGN are those with a mass larger 
than the Jeans mass in the conditions of the accreted matter and 
that have had sufficient time to collapse during the process. In 



T.J.-L. Courvoisier and M. Tiirler: Cascades of shocks in AGN 



7 



order to roughly estimate the collapse time of the clouds, we 
calculate the free fall of the cloud on itself 



tec — 



2GMj' 



where Sdoud is the initial size of the cloud while at large dis- 
tances to the active region and 



M.-(3)^( — ) 



(19) 



is the Jeans mass for a density of gas n at temperature T 
(see e.g. [Lang, Astrophysical formulae^ . The temperature of 
the matter far from the active region has been measured in the 
far infrared and found to cover a range with lowest measured 
temperatures around 50 K ( PoUetta et al. 2000l) . which we use 
here as representative of the accreted material. The density of 
the accreted matter is given by the mass accretion rate M, itself 
given by the Eddington luminosity and the accretion velocity 
Vacc- In a spherical geometry one finds: 



M 



Awr^Vac 



(20) 



This is likely to be a lower limit, as the large scale accretion 
flow is unlikely to be spherically symmetrical. The Jeans mass 
is thus expected to be smaller than estimated from Eg. 1191 

Using Eddington luminosity, an accretion efficiency of 
10%, and a black hole mass of lO'' solar masses, we show the 
Jeans mass at various initial distances to the active region as 
a function of the parameter x in Fig- El In the same figure we 
show in grey the region of the parameter space in which the 
cloud collapse time tec is shorter than the accretion time. Figure 
|2lthus shows, that for likely initial distances of lOpc and more, 
clumps of 100 solar masses and more can form and collapse, 
provided that the ratio x of the accretion velocity to free fall 
velocity is of several percent or less. 

Departure from the spherical symmetry, assumed here for 
simplicity, increases the density of the accreted matter and 
therefore decreases the Jeans mass at a given distance. Any ini- 
tial dumpiness of the accreted flow also favors the formation 
of clouds. It is therefore very likely that matter reaching the in- 
ner regions of a quasar will indeed be very clumpy rather than 
smooth. 

Neglecting self gravity and the cloud structure, we may 
estimate the deformation of the clouds while in the ac- 
tive region in the typical time between collisions. Following 
ICourvoisier et al. (1996)| the number of clouds N in a region 
of size R is given by the rate of collisions r: 



N = 



1 



^Kepler (R) 



20- 



1 



R. 



(21) 



c,15 



for a region of size R ~ 100 x Rs ~ 3 x 10^^ cm (correspond- 



py 



ing to Mbh = 10 M0), a cloud size Rc = 10 cm, and r, 
collisions per year. 

With the time between two collisions given by ^ and the 
shear generated by the Kepler velocity difference across the 



cloud one obtains that the cloud elongation during the time 
elapsed between 2 collisions along its orbit is about 



ARa 

(18) -r7 



20 



1 



1 



R, 



(22) 



cl5 



Since we neglected the self gravity of the clouds, this is an 
upper limit to the distortion of the clouds. Equation ( I22t shows 
that, while the clouds may be significantly distorted, they are 
likely to keep their identity and remain isolated on their orbits 
between major collisions. 

8. Summary and conclusion 

The picture that has been drawn is one of matter being accreted 
in the form of clumps of matter falling in the deep potential 
well of a massive black hole. These clumps interact with one 
another at first in the form of optically thick shocks and then, 
following an adiabatic expansion, through several optically thin 
shocks. The optically thick shocks give rise to the optical-UV 
emission, while the optically thin shocks give rise to the X-ray 
emission. 

The approach taken here has several elements that are sim- 
ilar to what would be expected from AGN models based on 
supernovae ( jAretxaga et al. 19^7 |, the main difference being 
that the source of energy is gravitational energy in the poten- 
tial well of the central black hole rather than nuclear energy or 
gravitational binding energy resulting from the collapse of sin- 
gle evolved stars. This then allows the individual events to have 
properties like total energy, event rate, or event duration that all 
depend on their environment, typically the density of clumps 
in the surroundings of the black hole or their angular momen- 
tum distribution. This dependence is important to understand 
the time variability of AGN as a function of their average lumi- 
nosity, as discussed in Torricelli-Campioni et al. (2000)| or the 
wide variety of observed properties seen in the class of Seyfert 
galaxies and quasars. 

The model we have developed here allows us to understand 
many of the properties of AGN, namely: 

- the temperature of the blue bump, 

- the small lag observed between the UV and optical AGN 
light curves, 

- the luminosity and shape of the hard X-ray emission, 

- the correlations between the blue bump and X-ray light 
curves, 

- the origin, velocity, and density of the broad line emission 
clouds. 

Moreover, the model does not make reference to any unspeci- 
fied source of viscosity. 

AGNs have blue bumps that differ considerably in impor- 
tance relative to other components of their continuum energy 
distribution. In the picture proposed here, this is related to the 
luminosity of the optically thick shocks compared to that of 
the optically thin shocks. It is indeed expected that a fraction 
of the matter (e.g. that which is only weakly bound in the ac- 
creted clouds) is not shocked first in optically thick conditions. 
Objects for which this fraction is large have weak blue bumps. 
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0.001 



0.01 0.1 

X = V^ccAff 



Fig. 2. The Jeans mass as a function of 
the distance to the black hole and the 
ratio X of the accretion velocity to free 
fall. Also shown in grey is the region 
in which cloud collapse is faster than 
accretion. Matter accreted from within 
the grey region of parameter space will 
be clumpy when arriving in the vicinity 
of the black hole, matter accreted from 
the white section of parameter space 
may be smooth, provided that it is ho- 
mogeneous to start with. 



The temperature of the blue bump remains, however, similar, 
as it is given by the Kepler velocity. 

There is one element of AGN physics that does not follow 
from the phenomenology we have considered, the existence of 
structured relativistic jets. Their generation does rely on fur- 
ther elements such as the rotation of the central black hole and 
magnetic fields. 
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